Strong Coupling Constant from the Photon Structure Function 
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We extract the value of the strong coupling constant a s from a single-parameter pointlike fit to 
the photon structure function F£ at large x and Q 2 and from a first five-parameter full (pointlike 
and hadronic) fit to the complete F 2 7 data set taken at PETRA, TRISTAN, and LEP. In next-to- 
leading order and the MS renormalization and factorization schemes, we obtain a s (mz) = 0.1183 ± 
0.0050(exp.)+°-™2s(theor.) [pointlike] anda s (m z ) = 0.1198 ± 0.0028(exp.)±g;°°^ (theor.) [pointlike 
and hadronic]. We demonstrate that the data taken at LEP have reduced the experimental error 
by about a factor of two, so that a competitive determination of a s from is now possible. 



The theory of strong interactions, Quantum Chromo- 
dynamics (QCD), is one of the corner stones of the Stan- 
dard Model of elementary particle physics. The precise 
determination of its fundamental parameter, the strong 
coupling constant a s , bears important implications for 
the validity not only of QCD itself, but also of even 
more fundamental theories, since these have to contain 
the Standard Model as an effective field theory in the 
low-energy limit. A more fundamental theory, which 
might explain the size of the strong coupling constant, 
has yet to be established. Therefore, a s must currently 
be extracted from experiment. Among the large variety 
of processes that have been used to this end, the most 
precise values have been obtained in Z-boson- and r- 
decays at LEP, scaling violations in structure functions 
at HERA, and quarkonium decay branching fractions and 
lattice calculations of quarkonium mass splittings, lead- 
ing - together with other, less precise measurements - to 
a current world average of a s {mz) = 0.1172 ± 0.0020 at 
the mass of the Z-boson, mz = 91.1876 GeV [0. 

When the photon structure function was first dis- 
cussed in the context of QCD, a precise determination of 
a s quickly emerged as one of its most interesting appli- 
cations. Due to the pointlike coupling of the photon to 
quarks, the leading order (LO, 0(a/a s )) § and next-to- 
leading order (NLO, 0(a)) [§ contributions to F^(x, Q 2 ) 
are calculable in QCD perturbation theory, if the virtu- 
ally Q in the deep-inelastic electron-photon scattering 
process is significantly larger than the asymptotic scale 
parameter A. Unfortunately, this pointlike contribution 
exhibits a power singularity at small Bjorken-a; [Q, which 
becomes rapidly stronger in higher orders The sin- 
gularity can be regularized with a non-perturbativc || 
or transverse-momentum |?]] cut-off, but then the sensi- 
tivity to a s is reduced and a dependence on the unphys- 
ical cut-off is introduced ||. It is then necessary to fit 
both a s and the cut-off to experimental data. Alterna- 
tively, the singularity can be canceled order by order in 
perturbation theory by retaining a hadronic boundary 



condition at a low starting scale Qq ||.In this case it is 
necessary to fit a s and the hadronic input to experimen- 
tal data. However, the evolution of the hadronic input 
to the physical scale Q is still predicted by perturbative 
QCD through inhomogeneous evolution equations Jl0| , 
and the negligibility of the hadronic input can be tested 
a posteriori. Both methods have been applied in the past 
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[|llj or a s {m z ) = 0.1081^. This value contributed to 
the world average in the 1988 fl§, 1990 fl3j, and 1992 
fTij issues of the Review of Particle Properties, but was 
then abandoned on the grounds that there were "no new 
results and the data do not contribute significantly to 
the average" p5[ . Since then it has been widely believed 
0, [l7], [L8], |lg |3 that the sensitivity of F 2 7 to a s is 
small. 

In this Letter, we point out that over the last decade 
a wealth of new F^ data has been collected at the e + e~- 
colliders TRISTAN and LEP, which extends to high av- 
erage values of Q 2 , (Q 2 ) < 780 GeV 2 . We demonstrate 
that the new data improve the sensitivity of F% to a s 
significantly and that a single-parameter pointlike fit as 
well as a five-parameter full (pointlike and hadronic) fit to 
PETRA, TRISTAN, and LEP data yield results, which 
are not only consistent with the world average, but also 
have competitive experimental and theoretical errors. 

We work in a fixed flavor number scheme with three 
active quark flavors (u,d,s). It is well known pl[ that 
for current measurements of Ftf the available hadronic 
energy squared W 2 — Q 2 (l — x)/x is not much larger 
than the production threshold 4)71^ of the heavy quarks 
(h = c, b, t), so that mass effects can not be neglected and 
the massive, fixed order O(o) expression for the Bcthc- 
Heitler process "f*{Q 2 )"f — * hh ]22| should be used instead 
of the massless, factorized 0(a/a s ) expression. For a 
consistent NLO (0(a)) analysis, we do not include the 
known (^3), but numerically small, 0(aa s ) corrections to 
the Bethe-Heitler process and omit the 0(aa 2 ) contribu- 
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tions from the process 7*(Q 2 )g — > ftft. The heavy quark 
masses are not well constrained from measurements of 
F2 ■ We adopt a charm quark mass of m c = 1.5 ± 0.1 
GeV in good agreement with recent precise determina- 
tions from threshold production at e + e~-colliders [ pi[ . 
At large Q 2 , charm quarks contribute up to 40 % to 
in the whole x-range, while at small Q 2 they contribute at 
most 10% below x ~ 0.2. The contribution from bottom 
quarks is suppressed by a relative factor 1/16 from the 
different quark charges, while the threshold for top quark 
production lies at extremely small x = 10 _5 ...10~ 3 , so 
that these contributions are both numerically negligible. 

Since we wish to omit spurious higher order terms, 
which arise from the convolution of NLO contributions 
to the parton densities with the NLO Wilson coefficients 
and lead to instabilities at large x ]2lJ , we work in Mellin 
moment space, where the convolutions reduce to sim- 
ple products, the evolution can be done analytically and 
without any approximations, and spurious higher order 
terms can be consistently omitted. The resulting pre- 
diction for F2 is then converted back to x-space using a 
fast inverse Mellin transform with logarithmic mapping 
and fitted to experimental measurements with the multi- 
dimensional minimization algorithm MINUIT p5| . The 
quality of the fit is measured in terms of the % 2 value per 
degree of freedom, % 2 /DF, for all selected data points. 
We include in our analysis all published measurements 
at the high-energy e + e~-colliders PE- 
M TRISTAN @ H, 0, and LEP 

||, H© EH E§M H ■ If more than one set of statisti- 
cally overlapping data exists, the most recent publication 
is used. We exclude from our fit the data published by 
the TPC/27 Collaboration at PEP §9], Q, since sev- 
eral data points, mainly at low x, are inconsistent with 
measurements published by PLUTO |27)) L3 p|, and 
OPAL in the range 1.9 < Q 2 < 5.1 GeV 2 Data 
where the charm component has been subtracted are also 
discarded. Statistical uncertainties and correlations be- 
tween data points due to the experimental unfolding are 
taken into account as provided by the experiments, while 
systematic uncertainties are assumed to be uncorrelated. 
Due to this assumption the values of x 2 /DF are expected 
to be on average slightly less than unity. If asymmetric 
errors are given by the experiments, the data points are 
taken at the center of the full error interval. Most exper- 
iments have not corrected for the finite virtuality of the 
target photon P 2 . We neglect P 2 in this analysis, since 
usually P 2 «Q 2 . 

For our pointlike fit, we identify the starting scale 
Qo with the asymptotic scale parameter A, so that the 
hadronic input vanishes automatically and only a single 
parameter (A, or equivalently a s (mz)) has to be fitted. 
As discussed above, this is only justified at large x and 
Q 2 , where the residue of the pointlike singularity is ex- 
pected to be small. Therefore, we perform our single- 
parameter pointlike fit only to a subset of data points 



TABLE I: x 2 /DF and a B (mz) values obtained in LO and 
NLO in the MS and DIS 7 factorization schemes with a single- 
parameter fit of the pointlike photon structure function . 
Also shown are the results obtained without LEP data and 



of F% collected 



TRA (26 



with very high 


Q 2 data. 






Scheme 


X 2 /DF 


a a (m z ) 




LO 

MS 
DIS 7 


7.9/ 19 
9.1/ 19 
8.1/ 19 


0.1260 ± 0.0055(e X ; 
0.1183 ±0.0050(e X ; 
0.1195 ±0.0051(e X ; 


+0.0061 / th N 
-0.0055 
+0.0029 ( , h] 
-0.0028 \'' n ) 
+0.0031 ( th N 
-0.0028 y''^) 


w/o LEP 
high Q 2 


3.2/ 7 
11.9/ 8 


0.1244 ±0.0126(ex; 
0.1159 ±0.0125(e X ; 


+0.0033 f th \ 
-0.0032 \'' n J 
+0.0018 / fh x 
-0.0018 \' jU ) 



with x > 0.45 and Q 2 > 59 GeV 2 . Very similar results 
are obtained with the widely used values of Qo = 0.5. ..0.6 
GeV @ |l|, (if |§, while choosing Q = 1 GeV sig- 
nificantly increases the value of % 2 /DF; two-parameter 
pointlike fits of a s and Qo are driven to Qo — A. In the 
first three lines of Tab. [j] we list the % 2 /DF and a s (mz) 
values obtained in LO and NLO. The NLO fit is per- 
formed in two factorization schemes (MS and DIS 7 |pl|) 
with different treatment of the pointlike Wilson coeffi- 
cient in F2 , but the numerical variation is found to be 
small. The total values of % 2 /DF as well as those for 
the individual data sets (not shown) lie around unity 
or below, indicating that the pointlike photon structure 
function and the fitted values of a s (mz) describe the 
data sets well within their statistical and systematic un- 
certainties. The experimental errors are determined by 
varying a s (mz) until the total value of x 2 is increased by 
one unit. To estimate the theoretical error, we vary the 
charm quark mass as indicated above and follow the com- 
mon convention of varying the factorization and renor- 
malization scales by factors of two about their central 
value, the physical scale Q. We then add these three in- 
dividual errors in quadrature. The LO value of a s (mz) 
is consistent with the NLO value within the expected ac- 
curacy, 0(a 2 ), and the theoretical error is reduced from 
LO to NLO as expected. In the fourth line of Tab. |[ we 
list the result of a fit without the LEP data. The ex- 
perimental error is more than doubled, showing that the 
LEP data have considerably increased the sensitivity of 
F2 to a s at high x and Q 2 . When data at all values of 
x, but very high Q 2 (Q 2 > 284 GeV 2 ) are fitted, the cen- 
tral value of a s (mz) remains virtually unchanged (last 
line of Tab. [j]). At very high Q 2 , the theoretical error 
drops by a factor of two, whereas the experimental error 
increases. Measurements of at a future linear e + e~- 
or e7-collider like TESLA at very high values of Q 2 and 
with small experimental errors will therefore lead to even 
more precise determinations of a s . 

The goodness of our pointlike fit may also be judged 
from Fig. [I], where the fitted data points are shown as full 
circles, while those that have been omitted from the fit 
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The Pointlike Photon Structure Function at Large Q" 



OPAL Q ! =59 GeV 2 
NLODIS y 

NLO MS 

LO 



OPAL Q-=76.4 GeV- 




TABLE II: Q , x 2 /DF, and a s (m z ) values obtained in LO 
and NLO in the MS and DIS 7 factorization schemes with a 
five-parameter fit of the hadronic photon structure function 
F%. Also shown are the results obtained without LEP data. 
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FIG. 1: Single-parameter fits of the pointlike photon struc- 
ture function^ compared to PETRA Q, TRISTAN g ^J, 
and LEP g ||, ||, || data at large Q 2 . The data 
points marked by open circles have not been used in the fits. 
Also shown is the hadronic contribution from a five-parameter 
NLO fit of the full photon structure function in the DIS 7 
scheme. 



are shown as open circles, and where the statistical and 
systematic errors have been added in quadrature. The 
theoretical curves are perturbatively stable, i.e. LO and 
NLO fits differ only by small amounts. The choice of 
factorization scheme affects the region of very large x, 
but it has only a minor effect on the description of the 
data. Also shown in Fig. |l] is the hadronic contribution 
from a five-parameter NLO fit of the full photon structure 
function in the DIS 7 scheme. It clearly falls from small 
to large x and Q 2 and amounts to only a few percent in 
the region that has been used in the pointlike fit. 

For our full (pointlike and hadronic) fit, we start from 
the observations that F% is dominated by the u-quark 
density in the photon and is only sensitive to the com- 
bined density of d- and s-quarks, whose contribution is 
furthermore suppressed by the smaller d- and s-quark 
charges. The gluon contributes to in LO only through 
a rather weak coupling to the quark singlet density in the 
evolution equations. A consecutive fit of the w-quark, d- 
and s-quark, and gluon densities shows, that only the 
first is well constrained by data and that the fit does 
not improve, when more degrees of freedom are added. 
Therefore we do not impose a hadronic boundary condi- 



Scheme 


Qo/GeV x 2 /DF a s (m z ) 




LO 

MS 
DIS 7 


0.79 ±0.18 121/129 0.1475 ± 0.0074(ex 
0.83 ±0.09 118/129 0.1198 ± 0.0028(ex 
0.85 ±0.09 115/129 0.1216 ± 0.0028(ex 


+0.0141 ( , h) 
-0.0072 V 111 ^ 
+0.0034 /.,-, 
-0.0046 V U1 V 
+0.0033 ( , h) 
-0.0050 I ln ) 


w/o LEP 


0.46 ±0.10 37/ 38 0.1147 ± 0.0047(ex 


+0.0282 
-0.0033 I U1 V 



tion for the gluon and assume, that the hadronic fluctu- 
ations of the photon are insensitive to the quark charge, 
i.e. we identify the hadronic boundary conditions for u- 
quarks and d- and s-quarks at the starting scale Qo- To- 
gether with a s (mz) and Qo, we then fit the parameters 
N, a, and (3 of our ansatz /J d+s (x, Qq) — Nx a (l—x) 13 to 
the full data set described above. In the first three lines 
of Tab. H we list the Qo, x 2 /DF, and a s (mz) values ob- 
tained with this five-parameter fit in LO and NLO. The 
starting scale Qo is perturbatively stable and is found 
to be close to the masses of the light vector mesons p, 
to, and 4> in contrast to earlier claims that the pertur- 
bative evolution of F% sets in only at rather high values 
of Qo ~ 2 GeV jll]]. The individual and total values of 
X 2 /DF lie again around unity or below, so that the fitted 
full photon structure functions describe the full data set 
well within the experimental uncertainties. Note that the 
X 2 value for the four TPC/2 7 points at Q 2 = 2.8 GeV 2 , 
which have not been used in the fits, is 18.0 and thus very 
large. The gluon density, generated with fg(x, Qq) = 0, 
turns out to be in good agreement with recent HI dijet 
data |^| . The experimental errors on the values of Qo 
and a s {mz) reflect an increase in \ 2 by one unit, when 
all other fit parameters are kept fixed. Due to the larger 
number of data points in the full fit, the experimental 
error turns out much smaller than in the pointlike fit. 
When the full fit is performed without the LEP data (last 
line of Tab. |l|), the experimental error is almost doubled, 
i.e. the impact of the LEP data is again impressive. A fit 
to LEP data only leads to almost identical results as the 
full fit. The theoretical error in LO and without the LEP 
data gets a large asymmetric contribution from doubling 
the factorization scale, which is highly correlated with an 
increase in the fitted value of Qo and which is drastically 
reduced in the full NLO fit. Similar results as those listed 
111 Tab. H are obtained, when only it-quarks are assigned 
a hadronic boundary condition. 

In Fig. U we compare our results to the fitted F^ data 
in the region of low x and Q 2 . This region is clearly dom- 
inated by the hadronic contribution and by the impact 
of the LEP data. A fit without the LEP data results in a 
rise of F% at low x, which is much too steep. The fits are 
perturbatively stable and the data are described almost 
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The Full Photon Structure Function at Small Q" 
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FIG. 2: Five-parameter fits of the fullphoton structure func- 
tion, compared to data from PETRA TRISTAN jjM], 
and LEP ||, ||, 0, ||, ||] at small Q 2 . The data points 
marked by open circles refer to the second experiment and/or 
Q 2 value. Also shown are the hadronic and pointlike contri- 
butions to the NLO fit in the DIS-v scheme. 



a s (m z ) = 0.1198 ±0.0054 (1) 

in NLO and the MS scheme, where the larger theoretical 
error has been added to the experimental error in quadra- 
ture. While our total error is slightly larger than those 
obtained in Z-boson- and r-decays at LEP, it is compara- 
ble to the errors obtained in deep-inelastic scattering at 
HERA and heavy quarkonium decays. This encourages 
us to combine our result with the current world average 
of 0.1172 ± 0.0014 to a new world average 

a s (m z ) = 0.1175 ±0.0014, (2) 

where the errors are assumed to be uncorrelated. 

In conclusion, we have for the first time fitted the now 
final PETRA, TRISTAN, and LEP data on the photon 
structure function K7 in NLO of perturbative QCD. Wc 
have extracted the value of the strong coupling constant 
oi s (mz) with competitive experimental and theoretical 
errors from a single-parameter pointlike fit to data at 
large x and Q 2 and from a five-parameter full (point- 
like and hadronic) fit at all x and Q 2 . Our analysis 
proves that the available F% data contribute significantly 
to a precise determination of a s and that future measure- 
ments of F% at linear colliders will have a large impact. 

We thank G. Kramer for many valuable discussions 
and a careful reading of the manuscript. S. A. and M. K. 
are supported by the Deutsche Forschungsgemeinschaft 
through Grant No. KL 1266/1-2. 



equally well in the MS and DIS 7 scheme. 



Since the total error on a s (mz) is smaller in the full 
fit than in the pointlike fit due to the larger number of 
data points, we adopt as our final result 
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